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Abstract. 

We investigate the quantum noise properties of a weakly nonlinear Duffing resonator 
in the deep quantum regime, where only few quanta are excited. This regime is 
dominated by the appearance of coherent multiphoton resonances in the nonlinear 
response of the resonator to the modulation. We determine simple expressions for the 
photon noise spectrum and find that the multiphoton resonances also induces a multiple 
peak structure in that noise. When the corresponding multiphoton Rabi oscillations 
are undcrdamped, zero temperature quantum fluctuations determine comparable 
populations of all quasienergy states which belong to a resonant multiphoton doublet. 
Most interestingly, the quantum fluctuations probe the multiphoton transitions by 
inducing several peaks in the noise spectrum of the resonator observablcs. In particular, 
the noise of the photon number contains complete information about the multiphoton 
states and their stationary populations via pairs of nearly symmetric peaks at opposite 
frequencies. Their widths are determined by the damping of the Rabi oscillations and 
their heights are proportional to the stationary noncquilbrium populations. A finite 
detuning from a multiphoton resonance generates a quasielastic noise peak at zero 
frequency. In addition, we relate the stationary populations of the quasienergy states 
with an effective quantum temperature and discuss the role of a finite temperature. 



PACS numbers: 03.65.Yz, 78.47.-p, 42.50.Hz, 85.25.Cp 
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1. Introduction 

Coupling a driven quantum mechanical oscillator to environmental fluctuations allows 
the oscillator dynamics to reach a stationary state. In the stationary state, energy 
is coherently absorbed from the pump and leaks into the environment via random 
dissipative transitions, which inevitably induce noise in the resonator. This occurs 
even at zero temperature where only environmental zero-point fluctuations (quantum 
noise) exist. The noise properties of a nonlinear oscillator determine many fundamental 
nonequilbrium phenomena such as quantum heating [H El E] and quantum activation 

i- 

Nonlinear oscillators are used as basic elements for quantum state detection or 
amplification. Examples of those are the Josephson bifurcation amplifier [5| IB] \7\ 18"! 15] [TU] 
and the cavity bifurcation amplifier [H [TTj . In this context, the noise properties of 
the resonator, which is used as detector or amplifier, determine the backaction of the 
measurement or amplification on the system itself [121 H31 E] • Clearly, it is desirable to 
keep the backaction as weak as possible, while on the other hand, a significant coupling 
of the amplification or measurement device to the system is useful in order to achieve 
a sufficiently strong detection or amplification efficiency. A fundamental lower limit for 
the introduced disturbance, however, will be set by the quantum noise. Hence, in order 
to design useful concepts for quantum state detection based on nonlinear resonators in 
the deep quantum regime, their quantum noise properties have to be addressed. 

The Josephson bifurcation amplifier takes advantage of the dynamically induced 
bistability due to the nonlinearity of the resonator. The eigenstates of the qubit are 
mapped onto the coexisting stable vibrational states of the resonator, which have 
different amplitudes and phase relations relative to the phase set by the external drive. 
Hence, they allow for a large discrimination power. Up to present, these amplifying 
devices mostly operate in a regime where many quanta in the resonator are excited. This 
implies that pure quantum fluctuations are typically small on average. Nevertheless, 
some experiments have been realized at low temperature where the relevant fluctuations 
are quantum mechanical in nature [H [10]. The regime of weak fluctuations has been 
the subject of intense theoretical investigation [U El El HH HH1 H51 HZ1 HBl HH]- It 
has been shown, that the stationary distribution over the quasienergy states of the 
driven oscillator at zero temperature has the form of an effective Boltzman distribution, 
thereby allowing to introduce the concept of an effective quantum temperature implying 
quantum heating even at T = 0. Signatures of the onset of quantum fluctuations 
can be seen in the relative intensities of the lines of the resonator noise spectrum 
[31 H5l [El [T71 [18] and in the appeareance of a fine structure in the spectral lines of 
resonators with comparatively large nonlinearities and large quality factors [31 EH]- The 
spectral fine structure yields detailed information on the quasienergy distribution [31 E] . 

Recently it has been proposed that nonlinear quantum detectors which operate 
in the regime of few quanta (deep quantum regime) would bring different advantages, 
such as a small back action, a large discrimination power with an enhanced readout 
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fidelity, and a sufficiently large measurement efficiency [19]. In the deep quantum regime, 
the frequency-resolved nonlinear response of the oscillator to the external driving with 
frequency cu cx shows a rich fine structure [U El [20], [2TJ [22] which is mainly generated 
by few-photon transitions in the resonator. The splitting of the typical Lorentzian 
resonance of a harmonic oscillator into a series of non-Lorentzian resonances and 
antiresonances reflects the intrinsic nonequidistance of the energy levels E n of a nonlinear 
oscillator. iV-photon transitions with the resonance condition — E = Nhu cx , 
N = 1,2, . . ., and the subsequent drift down along the ladder of the few-photon Fock 
states generate a pronounced nonequilbrium quasienergy distribution which is strongly 
different from the Boltzman-type [1] [2] [21]. Peaks or dips in the nonlinear response 
are a direct consequence of the nonequilbrium distribution over states with different 
oscillation amplitude and phase [1] [2] [2Tj . The signatures of such a characteristic non- 
Lorentzian lineshape of the response has been observed experimentally in a circuit-cavity 
QED set-up [23]. 

In this work, we investigate the noise properties of modulated nonlinear oscillators 
in the deep quantum regime. We consider the simplest example of a monostable 
anharmonic oscillator which has a quartic nonlinearity (Duffing oscillator). Such a 
weakly nonlinear Duffing oscillator has a remarkable symmetry: its energy levels E n 
with n < N are pairwise resonant for the same driving frequency w ex , E^_ n — E n — 
(N — 2n)hw cx . An example of the energy spectrum for the case N = 3 is sketched in 
Fig. [T] (a). After preparing the oscillator in its n-th excited state n < N, it displays 
periodic quantum oscillations between the n-th and the N — n-th excited states. During 
these oscillations, \N — 2n| photons are being exchanged between the oscillator and the 
modulation field. The oscillations of the photon number n are usually referred to as 
multiphoton Rabi oscillations. Their characteristic frequency, the Rabi frequency fi n jy, 
depends on the intensity of the driving field and on the number of photons exchanged. 
The Rabi frequency Qqn for the iV-photon oscillations is the smallest Rabi frequency. 
The multiphoton Rabi oscillations with N—n photons involved are underdamped if their 
Rabi frequency Q u n exceeds the dissipative rate of photon leaking into the environment. 
The latter is the oscillator relaxation rate 7. For 7 <C Qqn all the Rabi oscillations are 
in general underdamped. The periodically driven resonator reaches its stationary state 
on the timescale 7 _1 . 

In the stationary state, quantum noise induces - even at zero temperature - 
fluctuations in the photon number n. The dynamics of these fluctuations is characterized 
by multiphoton oscillations which manifest themselves as peaks in the noise spectrum 
S(cu) of n, located at plus/minus the Rabi frequencies tt n N- In the underdamped regime, 
the dissipative dynamics of the driven oscillator is most appropriately described in terms 
of random transitions between the oscillator quasienergy states. When the driving is 
resonant, the pairs of oscillator Fock states with n- and N — n-photons are resonantly 
superposed. The corresponding oscillator quasienergy states are a symmetric and an 
antisymmetric superposition of the two Fock states. Their splitting in quasienergy is 
given by the Rabi frequency fl n N- The corresponding peak in the noise spectrum at 
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(— )Q n N is due to random transitions from the state with (highest) lowest to that with 
the (lowest) highest quasienergy of the doublet. The peak intensity is proportional to the 
stationary occupation probability of the initial quasienergy state. Therefore, the noise 
spectrum offers a convenient way to directly probe the stationary distribution over all 
the quasienergy states. Moreover, for weak driving and exactly zero detuning from the 
multiphoton resonance, the noise spectrum of the n-photon transition is symmetric, 
i.e., S(oj) = S(—lu) and two inelastic peaks are signatures of an oscillatory decay 
of the fluctuations towards the stationary state. States belonging to a multiphoton 
doublet then have the same stationary occupation probabilities. For a weakly detuned 
modulation or a stronger driving, the spectrum becomes asymmetric. Besides, an 
additional quasielastic peak appears at zero frequency which represents incoherent 
relaxation of the fluctuations towards the stationary state. These features have some 
analogy in the spectral correlation function of a (static) quantum mechanical two-level 
system weakly coupled to a dissipative harmonic bath [23]. There, the spin correlation 
function is a sum of three Lorentzians. The two inelastic peaks are symmetrically located 
at finite frequencies and their width determines the inverse of the dephasing time. In 
addition, the quasielastic peak at zero frequency represents incoherent relaxation with 
the inverse relaxation time given by its width. In the driven system, the appeareance 
of a quasielastic peak depends on the intriguing interplay between the nonlinearity, the 
driving strength and the dissipation strength. 



2. Multiphoton Rabi oscillations of the Duffing oscillator 

We consider a periodically modulated quantum oscillator with mass m, eigenfrequency 
uq and a quartic (Kerr) nonlinearity described by the Hamiltonian 

H(t) = t^P 2 + 7^ muj oX 2 + ^a^ 4 + Fxcos(u cx t). (1) 

The modulation amplitude F is assumed to be so small that it induces only weakly 
nonlinear vibrations. This is guaranteed by the condition a A 2 <C mu^, with A(F) 
being the typical amplitude of the nonlinear vibrations. The modulation frequency u cx 
is chosen to be close to the oscillator eigenfrequency uo such that the detuning Su is 
small, i.e., 

5uj <C ojq , 5oj = ojq — u ex . (2) 

Our theory applies to hard as well as to soft nonlinearities a ^ 0, but for concreteness 
we will focus on the case of a hard nonlinearity, a > 0. 

The quantum dynamics of the weakly detuned and weakly nonlinear driven 
oscillator is most conveniently described in terms of the oscillator ladder operators a 
and o^, in a rotating frame determined by the unitary transformation 

R(t) = exp[— iu ex a i at] . (3) 

In the rotating frame, the typical time scale of the resonator dynamics is given by Su^ 1 , 
so that terms oscillating with frequencies ±2u ex and ±4cj cx average out and can be 
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Rotating frame 



(c) Quasienergy surface 
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Figure 1. Multiphoton Rabi transitions around the N = 3-photon resonance 
8u> = Suj^. In panel (a) we depict a sketch of the driving-induced resonant 3-photon 
transitions (red arrows) in the nonlinear oscillator. Likewise, the blue arrows indicate 
the noise- induced relaxation process. In panel (b), we show the underlying quasienergy 
spectrum as a function of the external frequency together with two zooms to the 
avoided crossings for the 2— and 3— photon Rabi transitions for v = 10~ 3 wo and 
/ = iy/10. In panel (c), we schematically indicate of the coherent multiphoton Rabi 
transitions (red arrow) and the dissipative transitions (blue arrows) on the quasienergy 
surface which results from a semiclassical approach, see text. The upper figure shows 
the less tilted case when / = v/10, while for the lower figure / = v. We emphasize 
that relaxational transitions at zero temperature typically occur in both directions, i.e., 
downwards and upwards along the quasienergy surface, which is in striking contrast to 
dissipative transitions in static potential surfaces, where only "downward relaxation" 
is possible. An escape due to "upward relaxation" is known as quantum activation [4]. 



neglected in the transformed Hamiltonian R(t) H(i) R'(t) — i HR(t) R^(t). Thereby, we 
obtain the RWA Hamiltonian 

H = 5un + vn(h + l)/2 + /(a 1 + a)/2, (4) 

where h = a^a is the photon number operator, v and / are the frequencies associated 
with the Kerr nonlinearity and the external field amplitude at the quantum scale 
xzpf = \/^/mwo, i.e., u = 3ax| PF /4^i and / = Fxzpf/V^H. In order to keep the 
notation compact we have set H — 1 in Eq. and in the remainder of the paper. 
The oscillator quasienergies e n and quasienergy states \ip n ) are the eigenvalues and 
eigenvectors of the rotating wave Hamiltonian, H\ip n ) = e n \ip n ). For vanishing driving, 
the quasienergy spectrum is given by 

s n = 5un + vn(n + l)/2 for / -»■ . (5) 
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We are primarily interested in studying the noise spectrum in presence of multiple 
multiphoton resonances E^_ n — E n = (N — 2n)co ex for n < N, or equivalently e^-n — £n 
for / — > 0. From Eq. ([5]) we find the resonant condition 

5lo = -v(N + l)/2 = 5u N . (6) 

Up to leading order in the driving, the quasienergy eigenstates \ip n ) for n < N ^ N/2 
are a resonant superposition of the pair of oscillator Fock states \n) and \N — n), i.e., 

\i, n )^{\n)±\N -n))/V2. (7) 

We choose the signs — and + for n < N/2 and N/2 < n < N, respectively. In the 
following, we refer to the resonant superposition of a pair of Fock states as resonant 
states or multiphoton states. The states which are not involved in a multiphoton 
transition (n > N and n = N/2 for N even) can be approximated as the corresponding 
Fock states \ip n ) ~ \n). The Rabi frequency VL n ^_ n of the multiphoton oscillations within 
the pair of Fock states \n) and \N — n) is given by the splitting of the corresponding 
levels e n and e^-n GSJ EI] as 



N-2n-l /Ar mi/2 



n ~\e e \-f( f \ {N ~ n) - 
U nN . n -\e n e N _ n \-J^-j ________ ^ 

The resonance condition in Eq. ([6]) is not renormalized by a finite driving within the 
RWA. Only for a comparatively larger driving / ~ ^ < wo, the multiphoton transitions 
have to be reinterpreted as tunneling transitions between semiclassical states [26j 127] . 

As we shall detail in Section HJ the multiphoton Rabi oscillations induce peaks 
in the spectral densities of oscillator observables only when the Rabi frequency Qqn 
for the multiphoton transition from the zero-photon ground state is larger than the 
noise-indcued level broadening of the relevant quasienergy levels Eq and En- In the next 
section, we will pave the wave for the calculation of the noise spectrum in this regime, 
by formulating the master equation for a weakly nonlinear oscillator and by evaluating 
the stationary occupation populations over the quasienergy states. 



3. Stationary dynamics in the deep quantum regime 

In the presence of a weak bilinear coupling to the fluctuations of a bosonic bath, the 
assumptions of small detuning and weak nonlinearity that underly the RWA naturally 
lead to a Liouville-von Neuman quantum master equation in the Lindblad form for the 
density matrix p of the weakly damped oscillator in the rotating frame [2T| [22] , 

'p = C P = -i[H,p\+ 1 {n + l)V[a}p + 1 nV[a ] }p 

V[0]p =([Op,0*\ + [0,pO*\)/2. (9) 

Here, C and T>[0] are the Liouville and the Lindblad superoperators, respectively. 
Moreover, 7 is the oscillator damping rate for which we assume that 7 C Wo. It 
results from a standard Ohmic bath spectral density J{oo) = 'yoj- In addition, n is 
the bosonic occupation number at frequency Uq and temperature T and is given by 
n = ( e w °/ fcsT - I)" 1 . 
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3.1. The stationary distribution 

For long times, the density matrix in the rotating frame p relaxes to a stationary state 
p°°, satisfying 

Cp°° = 0. (10) 

When the oscillator decay rate 7 is larger than the driving, 7 > /, the width of 
the resonant quasienergy levels e n induced by the bath fluctuations are larger than 
the corresponding Rabi frequency f2 n yv_„ of the multiphoton transitions. Then, the 
multiphoton resonances are smeared out and the coherent effects associated with 
multiphoton oscillations are strongly suppressed. Hence, dissipation sets a lower limit for 
the driving strength, / ^> 7, which has to be overcome in order to observe multiphoton 
transitions. On the other hand, for comparatively larger driving / ~ v, quantum 
fluctuations are significantly reduced and the oscillator is latched to a classical attractor 
at asymptotic times which are much larger than the typical relaxational time scale, 7" 1 . 
In this regime, the (quasi) stationary distribution of the population over the quasienergy 
states assumes Boltzman form [H El [31 [T71 [18] . 

Here, we restrict our analysis to the deep quantum regime where the driving is larger 
than the damping but smaller than the nonlinearity, 7 <C / <C v. Thereby, we have 
implicitly assumed a comparatively large nonlinearity v ^> 7. 



3.1.1. Underdamped regime: We start our discussion assuming that all Rabi oscillations 
are underdamped. Put differently, we assume that the smallest Rabi frequency Q n 
is larger than the relevant level broadening. We refer to this regime as the fully 
underdamped regime. Then, the off-diagonal matrix elements of p°° projected onto 
the quasienergy basis are negligible and we can set them to zero, i.e., we perform a 
secular approximation, 

p^ = (tpi\pP°M = for l^k. (11) 

Then, a balance equation for the stationary occupation probabilities p^ follows from 
Eqs. (jUJ) and (fTUl) according to 

-yiPu-^2^i,kPkk = o. (12) 

l^k 

Here, W^j- is the transition rate from state to state \fpi), 

Wi, k = 7 [(1 + n)|(^|a|V fc )| 2 + nKV^lVfc)! 2 ] (13) 

and 7/ is the width of quasienergy level Si given by 7; = Ylik^i Wi,k- We can now formulate 
more precisely the condition for underdamped Rabi oscillations to occur within the pair 
forming the narrowest resonance, which is VLqn ^> 70. 

The solution for stationary occupation probabilities up to leading order in the small 
parameters f/v and n is given in Ref. [21]: The pair of multiphoton states \ip n ) and 
\ipN-n) i n Eq. Q have equal stationary population, i.e., p^ n = p^_ nN _ n . The pair with 
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the narrowest resonance has the occupation probabilities = p^ N - The occupation 
probability grows algebraically with n < N/2 as 

N — n 

pZln+l = ^YPnn ^1 U < N/2 . (14) 

The states \ipi) with I > N have vanishing occupation probability, p^ = 0. As follows 
from the discussion above, the degeneracy pj^j = p™ N is approximate and is lifted for 
higher order in f/u. 

3.1.2. Quasienergy distribution close to a multiphoton resonance: One can easily 

generalize the above expressions to the case where the detuning 5u does not exactly 

match the resonant condition, 5u ^ 5un- Since the Rabi frequencies for the different 

pairs of resonant transitions in Eq. ([H]) are exponentially different, we can choose 

\Soj— Sun\ <C Qin-i, so that all the pairs of Fock states \n) and \N—n) with 1 < n < N/2 

are still resonantly superposed, except for 

9 9 9 9 

lV»o> =cos-|0>-sin-|JV> , and |^> = sin -|0> + cos -\N) , (15) 

with 9 = tan~ 1 [f2 A r /A r (5u; — 6u>n)]. The corresponding solution for the stationary 
density matrix close to resonance is 



Pnn = Po°°o tan 4 °-, pTx = (&N tan 2 t , 
N — n 

p n ° +ln+l = —^-p™ n for l<n<N/2. (16) 

3.1.3. Partially underdamped regime: Next we consider a comparatively large 
relaxation rate 7, so that the narrowest Rabi resonance is overdamped but the remaining 
resonances are still underdamped, VLqn <C Nj <C fiuv-i. We refer to this regime as 
the partially underdamped regime. Then, incoherent multiphoton transitions from the 
ground state |0) to state \N) with a small rate Qq N /(Nj) and the subsequent emission 
of excitations into the bath determines a small but finite occupation of the resonant 
states p™ n , n > 1. Formally, the stationary distribution p°° can be obtained by setting 
all the off-diagonal elements of p^? to zero except for p^ and p^ N and solving Eq. fTTOj) . 
Thereby, we find 

(0|pP°|0>«l, (N\p°°W) « njW(JVV) 

p™ = f4v/(iV7 2 ) for fijWW) » exp[-u /(k B T)] 
N — n 

P Zm + i = ^YPnn for \<n<N/2. (17) 

The crossover between this solution and the fully underdamped solution Eq. (fl4|) is given 
in Ref. [21J. Both stationary nonequilbrium distributions are determined by quantum 
fluctuations and are very different from the equilibrium Boltzmann-type distribution 
when a driven resonator is latched to a classical attractor. 
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3.2. The nonlinear response of the oscillator 

In the steady state regime, t 3> 7 _1 , the oscillator state is described by the time- 
independent density matrix p°° in the rotating frame and the oscillator dynamics is 
embedded in the time-dependent reference frame R(t). The mean value of an observables 
O is 

(O(O)oo = lim (0(t)) = Tr[p°°tf(t)OR(t)} . (18) 

t— >oo 

Therefore, the stationary oscillations of the position expectation value {x(t))oo are 
sinusoidal, 

(x(t))oc = v^xzpf cos (w cx t + v9)|(a) 00 |, (a)oo = y~]pul{'ij>i\a\'ij)k) ■ ( 19 ) 

Ik 

It has been shown that the nonlinear response (a; of the oscillator as a function of 
shows resonances and antiresonances in the deep quantum regime [20| |2T| [22] . The 
response is proportional to the transmitted amplitude in a heterodyne measurement 
scheme and it has already been measured for a weakly nonlinear oscillator [23] . Clearly, 
such a measurement scheme, or more general, any measurment scheme which probes 
stationary mean values as opposed to correlations does not allow to resolve the different 
degenerate resonances separately. Neither, they allow us to access the stationary 
distribution directly. This becomes possible only when correlations, e.g., via noise 
spectra are measured. In the next section, we show that this can indeed be achieved by 
measuring the spectrum of the photon number noise. 

4. The noise spectrum in the deep quantum regime 

4-1. Definition of the noise spectrum 

The Lindblad master equation ([9]) in general also allows to investigate transient 
phenomena and correlation functions. Its formal solution for a given initial state po 
is given by pit) = e ct po- Moreover, a general correlator (O'(t')O(t)) can be evaluated as 
the mean value of the operator O' at time t' with the virtual operator R) {t)0 R{t) p{t) at 
time t. This view has been established several decades ago by the Lax formula [281 E2] 
according to 

(O'(t')O(t)) = Ti{W{t')0'R{t')e c ^^W{t)OR{t)p{t)} (20) 
For long times t 3> 7 _1 , we find that 

(0(t + St)C(t)) oo = lim (0(t + 5t)0'(t)) 

t— >oo 

= Tr{R\t + 5t)OR(t + 5t)e c5t R ] \t)0' 'R(t) p°°} . (21) 

In general, such correlators are periodic functions of the preparation time t. The noise 
spectrum is defined as a double average over quantum fluctuations and the time t. 



Quantum noise properties of multiphoton transitions in driven nonlinear resonators 10 



Here, we are specifically interested in the noise spectrum S(u) of the autocorrelator 
of the photon number n, (n(t + 5t)n(t))<x,. From Eqs. and (12~T1) . we find 

(n(t + tffc)n(t))oo = Tr{ne £<5 *np°°} . (22) 

Since this correlator does not depend on the initial time t as a consequence of the RWA, 
we can define the noise spectrum in terms of a single average over quantum fluctuations 
according to 

POO 

S(u) = 2Re dte^(n(t)n(0))oo. (23) 



JO 

It is useful to separate the contributions to S(u) into those coming from the expectation 
value of n, and those from its fluctuations, i.e., 

POO 

5S(u) ee 2 Re/ dte^iSntySniO))^. (24) 



JO 

Here, 5h is the operator for the photon number fluctuations, i.e., 5n = h — (n)^. 

Our path to compute the noise spectrum consists in three steps: i) We express the 
virtual preparation hp°° in terms of right eigenvectors of the superoperator C. ii) We 
plug the resulting decomposition into Eq. fl22|) . Then, each term decays exponentially 
with a different exponent which is given by the corresponding eigenalue of C iii) We 
compute the Fourier integral in Eq. (|23|) . which thereby yields a sum over (overlapping) 
Lorentzians. 

The general expression, which is useful for a concrete numerical evaluation, for 
the noise spectrum given in terms of the eigenvectors and the eigenvalues of £ is 



derived in|Appendix A[ In the next section, we consider the special case of underdamped 



multiphoton Rabi oscillations. 



4-2. Noise spectrum in the underdamped regime 

When all the multiphoton Rabi oscillations are underdamped, Q^o ^> T^, the 
coherences IV'iV-nKV'nl and \i>n){i>N-n\ are approximate eigenvectors of the Liouvillian 
C. Then, 

C\l/)n)(i>N-n\ = - (r n - iQ nN - n )\iJ n )(tp N - n \ for U < N/2, 
£|^jV-n)(^n| = - (^n + i^nN-n)\^N-n)(^n\ for U < N/2 . (25) 

with the level widths being given as T n = 7 n = 7(77, + 1/2)N + 777. for n < (N — l)/2. 
For N odd, F(jv_i)/2 = 7(1 + 2n)(5A r + l)/8 + 772. Up to leading order in f/u, the 
decomposition of the virtual preparation np°° in terms of right eigenvectors of £ has 
the simple expression 

np°°*(N/2)p°°- E (N/2-n)pZ(\^n){^N-n\ + \^N-n)M) ■ (26) 

n<N/2 

Clearly, each term of the above decomposition yields a Lorentzian peak in the noise 
spectrum S(u). The first term yields the contribution to S(cj) from the expectation 



Quantum noise properties of multiphoton transitions in driven nonlinear resonators 11 



value of n, (N/2) 2 5(u). The remaining terms yield inelastic peaks associated to random 
transitions between quasienergy states belonging to the same multiphoton doublet. 
Since the populations p^ n and Pn-uN-u are approximately equal, peaks at opposite 
frequency have approximately equal intensity. By putting together Eqs. ( 12"2"j) . ( 12"3"j) . 
(ESD, and (J26D, we find S(u) = (N/2) 2 5(uj) + 6S(u) with 

5S(u)k S n (u) + S N „ n (u), (27) 

n<N/2 

9 (u) - 9 ( u) - 2T "Pnn(N/2-n) 2 

S n (u) - S^i-u,) - {u _ nnN _ n? + n ■ (28) 

Hence, the Lorentzians are centered at the multiphoton Rabi frequencies fl n N- n and 
have a resonance width of T n . The factor (TV — 2n) 2 /4 is the leading order expression 
for the squared matrix element \{ip n \n\i^N-n)\ 2 ■ Remarkably, the line intensities depend 
only weakly on the driving / and on the temperature through the stationary distribution 
p™ n . Up to leading order, the driving / enters only in the splitting of the lines through 
the Rabi frequencies. Notice that Eq. (127)1 is valid only in the vicinity of a multiphoton 
peak since terms of order 7 are not taken into account. In order to evaluate the tails of 
the peaks more precisely, one has to take into account the contribution stemming from 



all eigenvectors of £, see Appendix A 



In the left and right panels of Fig. [21 we show the noise spectrum S(oj) for the 
cases N = 2 and A^ = 3, respectively. The noise spectrum for N = 2 shows a pair 
of symmetric peaks which correspond to the transitions j^o) "H- \^2}- Likewise, the 
noise spectrum for A^ = 3 displays two pairs of symmetric peaks corresponding to the 
transitions \ipo) -B- ^3) and l^)- The green dashed lines mark the results from 

our approximate analytical formula in Eq. ( |27j) while the yellow solid lines show the data 
obtained by numerically evaluating the expression in Eq. (IA.2j) . An excellent agreement 
is found. 

In Fig. 12k), additional smaller side peaks of the order of f jv are also visible, see 
the gray lines representing a ten-fold zoom. They are not associated to any resonant 
transition between multiphoton states and are thus not captured by the leading order 
expression given in Eq. ( I2"7j) . The particular subleading peaks in Fig. [2^) belong to the 
transitions |^o) 

These features have a direct analogy in the spectral correlation function of a static 
quantum mechanical two-level system which is weakly coupled to a dissipative harmonic 
bath [21]. For a general biased two-state system with anticrossing energy levels, the 
pair correlation function is a sum of three Lorentzians. The two inelastic peaks are 
symmetrically located at finite frequencies and their width determines the inverse of the 
dephasing time. For a biased static two-level system away from resonance, an additional 
quasielastic peak at zero frequency appears which represents incoherent relaxation with 
the inverse relaxation time given by its width. Since we consider here the case strictly 
at resonance (in the RWA), no zero-frequency peak is present. 
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Figure 2. Photon noise spectrum at the N = 2- and the N = 3-photon resonance 
are shown in panel (a) and (b) respectively for v = 10 _3 wo- / = 7 = f2o2/10 

in (a) and 7 = Q03/IO in (b). Shown are the approximate results obtained with Eq. 
(|27p (dashed green lines), and the results from a full numerical solution of the general 
expression for the spectrum derived in Appendix A (orange solid lines). The gray lines 
in panel (a) mark a zoom to the subleading off-resonant transitions. 



4-3. Photon anti-bunching 

In general, the photon emission characteristics of a quantum mechanical resonator 
can show peculiar nonclassical features. For instance, counterintuitive correlation 
phenomena such has photon antibunching can occur, where the photon number 
correlation function for short delay times is smaller than the one for classical, 
uncorrelated photons. This implies that the probability for photons to arrive in 
pairs is suppressed [30]. Our approach provides a natural framework to investigate a 
possible non-Poissonian statistics of the multiphoton events in the nonlinear resonator. 
Therefore, we consider the normalized photon number correlation function or second- 
order coherence function defined as 

(2) (at(t) a t(t + r)o(t + ^0(^)00 

9 1 ) (at(i)a(t)) 0O (at(t + r)a(t + r)) 0O - 1 ) 

For long delay times r, the counts of two photons with a delay time r are statistically 
independent events, g^(r —> 00) = 1. For vanishing delay times, we have 



gM{r^0) = 1+ v " /Q ° Z;; 00 (30) 



(n 2 )oo - (n)^ - (n) 
(n) 2 



1 00 



Photon antibunching corresponds to the case g^ 2 '{r = 0) < 1. For the fully 
underdamped case, we find the expression 

m , s 2iV(N-l)+4Vf- 1 n(n- jV)p nn 1 , . 

9 {2 \r = 0) = — ^ >- ^ ^ ^ = 1 - - , (31) 

which represents the known result of the second-order correlation function of the 
electromagnetic field [3U]. Hence, the oscillator displays photon antibunching close to 
a multiphoton transition. The second-order coherence of the stationary state of the 
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quantum Duffing oscillator at the iV-th multiphoton resonance has the same value as 
the second order coherence for an oscillator prepared in the single Fock state \N), in 
spite of its fluctuations over the quasienergy states. 



4-4- Lineshape of the noise spectrum close to a multiphoton resonance 

In presence of a small detuning from the multiphoton resonance, 5u — 5ujn ~ &on, 
the states \ipo) and \iPn) are no longer a resonant superposition of the Fock states |0) 
and \N). Hence, the corresponding stationary occupation probabilities p o an d Pnn, 
given in Eq. (Tl6|) . become significantly different. In turn, the pair of peaks So(u) and 
Sn(u), which are associated to the transitions |^ ) I^n), become asymmetric such 
that Sq(oj) 7^ Sn{— oj). This behavior is shown in shown in Fig. Hk) for the case around 
the 3-photon resonance. The peak lineshapes can readily been evaluated and we find 

2r p °giV 2 (sinflcosfl) 2 2r^ jV 7V 2 (singcosg) 2 
S0{UJ) ~ F^FW' N{U) " ^-eo + e^ + Tl - (32) 
Their distance increases with the quasienergy splitting, — £o — sgn(5u — 5un)(QI n + 
N 2 \§tu — Scun] 2 ) 1 ^ 2 , whereas the peak width does not change close to the multiphoton 
resonance, 5u — Sun ~ Qqn- The asymmetry is determined by the stationary occupation 
probabilities yogg and p^ N - From Eq. (IT61) . we find 



S{ U ) _ pSS _ coll <L 



S(-uj) p% N 



n 



ON 



\£n — £o| — N(5u — 5u 



N, 



1 



(33) 



The above expression is valid for u close to the center of the largest peak, u ~ En — £o, 
and \5cu — Sun\ not too large such that S(±u) 3> 7. 

In addition to the peaks at finite frequencies (which induce decaying coherent 
multiphoton Rabi oscillations), also a zero frequency peak appears. This quasielastic 
peak is associated to incoherent relaxational decay of the multiphoton Rabi oscillations 
and is also known for the noise correlation function of a static biased quantum two level 
system [23] . 

In Fig. |3b), we show the logarithm of the asymmetry ratio given in Eq. fl33|) . The 
asymmetry shows a clear maximum at approximately — Eq. 

To further illustrate the asymmetry in the peak heights, we show in Fig. |3b) the 
peak maxima associated to the transitions l^o) — * IV^) an d l^) —> l^o)- At the 3-photon 
resonance (black dashed vertical line), both peaks are equal in height (symmetric noise 
spectrum). Away from the resonance, the low (high) frequency branch aquires more 
spectral weight for negative (positive) detuning. 



4-5. Photon noise at zero frequency 

Fluctuations of an oscillator (quasi)energy induce a broad (with width oc 7) zero 
frequency peak in the noise spectrum of an observable whose mean value depends on the 
(quasi)energy [31]. For weak driving / <^v and at a resonance \5u — 5ojn\ <§C ^W, the 
quasienergy states of the Duffing oscillator have large fluctuations as several quasienergy 
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Figure 3. (a) Asymmetric structure of the photon noise spectrum at frequency 
5ui = 5ui3 + 5, i.e., out of resonance for a detuning S = 1.6 x 10 -4 z^ for the same 
parameters used in Fig. [2b) (orange solid line). In addition, we show in the background 
the symmetric photon noise at the resonant frequency Suj^ (grey shadowed area). 
Moreover, we depict the inverted case S — > —5, which shows a symmetric behavior 
under the reflection lu — > —u (green solid line), (b) Noise asymmetry via the logarithm 
of Eq. ((33|) for the same parameters as in (a), (c) Height of the photon noise peak for 
the transition |V>o) — > 1^3) (orange solid line), and — > |^o) (green solid line) as a 
function of the external frequency. The peak maximum is located at <5w3 ± S. 



states have comparable occupation probabilities even at T = 0. However, the mean 
value of n becomes independent from the quasienergy, (V'nl^l'i/'n) « N/2 for n < N . As 
a consequence, the contribution to the noise spectrum of h coming from fluctuations 
SS(oj) does not have a peak at zero frequency since 5S(0) oc 7. Close to resonance, 
when \5cu — 5un\ ~ &on, two dynamical effects compete: on one hand, the quasienergy 
fluctuations quickly decrease for increasing detuning, i.e., moving away from resonance 
as the occupation probability of the state \ipo) approaches one. On the other hand, the 
mean value of n becomes strongly dependent on the quasienergy. As a result of this 
competition, the intensity of the zero frequency noise plotted as a function of Sou has 
two maxima at the two opposite sides of the resonant value Sun- In Fig. HJ we show 
the zero frequency noise for the special case N = 2. The yellow solid line represents 
the intensity at zero frequency computed numerically, while the green dashed line is the 
leading order contribution (in f/v) 

ob\U) — 0) ~ 47 sin - cos - sin cos - 

v ; 1 2 2 V 2 2 




4.6. Noise spectrum towards the semiclassical regime 

Next, we investigate the noise spectrum for larger driving strengths, / < v. In order to 
illustrate how the noise spectrum changes for increasing driving, we show the intensities 
of the brightest peaks as a function of the driving strength for the N = 5-photon 
resonance in ; see Fig. [5^). In Fig. |5b), we also show the quasienergy spectrum, and the 
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Figure 4. Photon noise at zero probe frequency as a function of the external frequency 
for the parameter set used in Fig. [^a) evaluated around the second multiphoton 
resonance, 5lo ~ —3v/2. Shown is the comparison of the approach Eq. (|34[) as a green 
dashed line with the numerical simulation as a orange solid line. 



noise spectrum for a comparatively large value of the driving amplitude / = v is shown 
in Fig. |5b). A peak in the noise spectrum at frequency uj = E\ — e k is associated to a 
single transition — > and is given by 



Hence, the relative intensities of a pair of peaks at opposite frequencies is still related to 
the occupation probability of the corresponding initial states through S(ei — e k )/S(e k — 
£1) = Pkk/PU- 

For weak driving, we have three pairs of approximately symmetric peaks as 
described by Eq. (I27p . Each peak corresponds to a transition between two states 
belonging to a multiphoton doublet of quasidegenerate states: \ipo) -H- \ipi) \1jj4), 
and ^2) ^ 1^3)- F° r increasing driving, the spectrum becomes increasingly asymmetric. 
For moderate values of the driving, the noise spectrum undergoes two major qualitative 
changes: i) the peak at zero frequency becomes clearly visible; ii) a pair of peaks 
corresponding to the transitions |?/>i) -H- IV'3) acquires a significant intensity. For / = v, 
the peak associated with the transition ^3) — > is even the second brightest peak. 

These qualitative changes can be explained in terms of a semiclassical description 
valid beyond the weak driving limit. The RWA Hamiltonian in Eq. (j4j) can be rewritten 
in terms of rotating quadratures, and interpreted as a quasienergy surface in phase space 
[26. 27 . It lias the shape of a tilted mexican hat and is sketched in Fig. Q] (c) for two 
values of /. The larger / is, the stronger is the induced tilt. The local maximum and 
the minimum of the quasienergy surface are the classical attractors. In the static frame, 
they describe stationary oscillations with a small and a large amplitude, respectively In 
the vicinity of the attractors the vibration amplitude and the slow part of the oscillation 
phase display slow vibrations with frequency oc 5u. In absence of resonant transitions, 




Ik 



2 P ™\{^ l \rfa\^ k )\ 2 ( 1 (a l - a k ) 2 + Y l + r fc ) 
(w + ex - Sk) 2 + (l(ai - a k ) 2 + r, + T fc ) 2 ' 



(35) 
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Figure 5. (a) Height of the photon noise peaks corresponding to the transitions 
withiin the pair \ip n ) <-> \ip m ) as a function of the driving strength /. Each pair is 
marked by a solid and a dashed line in the same color. In addition, we depict the rise 
of the zero frequency peak (black solid line) as the driving strength increases. The 
black horizontal lines indicate the expected values of the noise level evaluated up to 
leading order in f /v by using Eq. ([27]) . The parameters are v = 10~ 3 wo, Suj = SLO5. In 
panel (b) , we show for the same parameters the quasienergy spectrum as a function of 
the driving strength /. In panel (c), the photon noise spectrum as a function of the 
probe frequency u is shown for a large driving strength / = v. 



each quasienergy state can be associated to a quantized quasiclassical orbit which lies 
on the internal surface around the local maximum, on the external surface, or along 
the quasienergy well around the minimum. For very weak driving, / <^ vj \f2{N + 1), 
the quantum mechanical Fock states \n) with n < N/2 are associated to quasiclassical 
trajectories on the internal surface around the local maximum, whereas the Fock states 
with photon number n larger than N/2 are associated to semiclassical orbits on the 
external surface. Within this representation, the multiphoton transitions can then be 
reinterpreted as tunneling transitions between the internal and the external parts of the 
surface [261 EZ] • For comparatively larger driving, the zero-point quasienergy associated 
to the slow vibrations around the minimum (oc 80S) becomes smaller than the dynamical 
barrier height. Then, quasienergy states appear which are localized in the quasienergy 
well. In turn, the noise spectrum becomes qualitatively different from the one for weak 
driving. The small quantum fluctuations around the minimum of the quasienergy surface 
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can be described in terms of an effective auxiliary oscillator with ladder operators b and 
b' and are given by 

a = ah + b cosh — ¥ sinh r* h . (36) 

Here, is the amplitude of the stationary oscillations rescaled by v^zpf [T7J IT8] . 
They can be mimicked by a local effective quantum temperature T e = (2kg lncothr^)" 1 
which depends on the squeezing factor r* h [H El [TTJ, [18]. For / = u, the states \ip2), 
|^ 3 ), and \ipi) can be identified with the groundstate and first two excited states of the 
auxiliary oscillator (but in the remainder of this discussion we keep the same labels for 
the states as in the weak driving limit). The level spacing e 3 — e 2 is of the order of the 
frequency of the slow classical oscillations of the amplitude and slow part of the phase. 

Such oscillations appear in the noise spectral density of a classical oscillator as a pair 
of peaks. In a nonlinear quantum oscillator whose quasienergy levels are not equidistant 
and their distance exceeds the damping strength, the classical peaks have a "quantum" 
fine structure [3J. In the present case of the Duffing oscillator, the classical noise peak 
is splitted into two peaks associated to the nearest neighbor transitions between the 
ground state and the first excited state, and the first and the second excited state, 
\ip2) 1^3) an d \4>z) l^i), respectively. Their peak height is proportional to the 
square of the rescaled vibration amplitude ah and to the occupation of the initial state 
p™ n . The latter, in particular, is governed by the quantum temperature T e . For the ratio 
of the peak heights, we find [3] 

~ ^ ~P%~ coth 2 r * „ P§§ „ g(£l ~ gg) (37) 

Next nearest neighbor transitions can also yield peaks in the noise spectra of a Duffing 
oscillator [18]. In the present case, the transitions ^2) yield a pair of dimmer 

peaks, however, located at frequencies outside the frequency range shown in Fig. 

In the weak damping, weak driving regime discussed so far, the quasienergy well 
around the minimum is still very shallow, and the oscillator can escape from the small 
amplitude attractor via tunneling. Therefore, the oscillator is not latched to any of the 
attractors and the noise spectral density has also peaks which are associated to intrawell 
transitions. In particular, the pair of peaks with the smallest splitting describes coherent 
tunneling oscillations between the internal and the external part of the quasienergy 
surface (coherent dynamical tunneling or multiphoton Rabi oscillations). 

Before closing this section, we mention that for the stronger driving / = v, also a 
zero frequency peak appears in the noise spectrum, see Fig. [Hb), although the frequency 
detuning has been fixed to the 5-photon resonance 5u = 5u^. However, as discussed 
above, this resonance condition is only valid for small / < which is obviously not 
fullfilled. So the larger driving induces an effective small detuning away from the exact 
avoided quasienergy level crossing and generates an effective bias. Then, a relaxation 
pole appears in the relevant self energy [21] which corresponds to a quasielastic relaxation 
peak at zero frequency. 
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Figure 6. Photon noise at the third multiphoton resonance, Suj = Stu^, as a function 
of the probe frequency. In panel (a), we show the behavior when going from the 
coherent to the incoherent regime by increasing the damping constant from 70 <C ^03 
to 73 = fio3 ( a t T = 0). In panel (b), we fix the damping constant to 7 = f2o3/10 and 
show the noise temperatures increasing from zero (T = 0) up to finite temperatures 
(T ^> u>q). The remaining parameters are v = 10 _3 a;o and / = v/10. 



4-7. Dependence of the noise spectrum on damping and temperature 

So far, we have analyzed the case of zero temperature and small damping, n < 1 and 
7 <C ^otv- In this section, we briefly address how the noise spectrum is modified for 
larger damping and finite temperature by presenting numerical results of the spectrum 
in a broad parameter range. 

In Fig. we show S(ou) for different values of the damping for the 3— photon 
resonance where 5u = 5u^. As expected, the peaks in the noise spectrum get broader 
for increasing damping. Outside the fully underdamped regime, the two peaks of the 
pair associated with the transitions |^o) ^ IV^) start to overlap and eventually merge 
into a single peak at zero frequency. Thereby, the zero frequency noise is no longer 
suppressed S(u ~ 0) oc 7 -1 , since incoherent relaxation prevails over coherent decay for 
large damping. The peaks associated with the underdamped transitions \ipi) <H- \1jj2) 
are still described by Eq. ( 127|) . even when the spectrum has a peak at zero frequency. 
The decrease in the peak intensities reflects the decrease of the populations and 
in the partially underdamped regime. 

The dependence of the noise spectrum on temperature is shown in Fig. [Hb) and 
behaves qualitatively similarly. For small temperatures n < 1, the spectrum is described 
by Eq. ( )27|) . The temperature dependence enters in the line widths of the quasienergy 
levels as well as in the stationary distribution p™ n . For larger temperatures, the two low- 
frequency peaks merge into a single peak at zero frequency and the side peaks becomes 
increasing broader as expected. 
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5. Conclusions 

In recent years, the rich phenomenology of driven and damped nonlinear quantum 
oscillators has been impressively consolidated, including their nonlinear response 
behavior in form of resonant and antiresonant amplification, quantum coherent 
multiphoton Rabi oscillations, quantum activation and quantum heating. Gradually, 
the nontrivial effects visible in noise correlation functions have also moved to the focus 
of interest. Those become relevant whenever a nonlinear quantum oscillator is used 
as a central element in an amplifier or quantum measurement device. In this work, 
we have analyzed the noise properties of the quantum Duffing oscillator in the regine 
when only few quanta are excited. Then, the nonlinear response shows pronounced 
multiphoton peaks which are associated to resonant multiphoton Rabi oscillations. 
The noise properties of these multiphoton transitions show a rich phenomenology. To 
obtain the noise spectrum by analytical means, we invoke the Lax formula for the 
autocorrelation function of the photon number at different times and calculate its 
Fourier transform. Exactly at a multiphoton resonance, the noise spectrum consists 
in a collection of pairs of related resonances which are located at opposite frequencies 
and which are equal in height. Each pair is associated to a multiphoton doublet. In 
spite of large fluctuations over the oscillator quasienergy, no quasielastic peak occurs at 
zero frequency. This is a consequence of a special symmetry of the quantum Duffing 
oscillator: all quasienergy states which are associated to a multiphoton doublet have 
the same mean value of the photon number n. 

Slightly away from a multiphoton resonance, the noise spectrum becomes 
asymmetric and the two resonances are no longer equal in height. In addition, as the 
mean values of n become different for quasienergy states with comparable occupations, 
the quasielastic peak emerges. Since the quasienergy fluctuations are suppressed away 
from a multiphoton resonance, the intensity of the quasielastic peak as a function of the 
detuning displays a maximum at the two opposite sides of the resonant value 5un- 

Our approach also allows us to evaluate the transition to the semiclassical regime 
by increasing the photon number by a larger driving amplitude. Then, a quasiclassical 
quasipotential landscape in phase space is a convenient tool to understand the stationary 
nonequilibrium dynamics. This view directly leads to quantum mechanical squeezed 
states which exist close to the local minimum of the quasienergy landscape. A harmonic 
expansion allows us to characterize the quantum fluctuations via an effective quantum 
temperature. At larger (real) temperature and damping strengths, all these quantum 
coherent features are washed out. 

Although the time-resolved detection of noise properties of quantum observables 
of driven resonators requires considerably more experimental effort, we are confident 
that future experiments will soon elucidate the importance of quantum noise in these 
systems. 
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Appendix A. Numerical evaluation of the noise spectrum 

In order to numerically compute the noise spectrum S(u), we diagonalize the Liouvillian 
superoperator C, 

Ck R = X k k R , k L C = X k k L (A.l) 

Here, {/cl} and {k R } are the sets of the left and right eigenvectors of C, respectively. 
They constitute a biorthogonal system implying that trjA^A;^} = S kk i. The 
corresponding eigenvalues are in general complex. The stationary density matrix p°° 
is the only eigenvector with zero eigenvalue. All other eigenvectors represents transient 
dissipative processes. The corresponding eigenvalues have a negative real part. Note 
that we have implicitly assumed that £ admits a spectral decomposition. Put differently, 
we have ruled out that any Jordan block has dimension larger than one. 
With this at hand, the noise spectrum in Eqs. (123]) can be rewritten in terms of fcj,, k R 
and X k as a sum of partial spectra characterized by Lorentzians according to 

at. ,\ _ 2E fc (tr{nfc R }tr{fc[np 00 })ReA fc 

The sum extends over those eigenvectors k R which belong to non-zero eigenvalues X k . 
Thereby, we have not included the elastic Rayleigh peak (n) 2 ^^) which trivially comes 
from the stationary state k R = p c 
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